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ABSTRACT 

The problem of thermalization of a quark - gluon plasma is addressed in the frame- 
work of thermal field theory. Within a simple approximation, the full quantum 
relaxation problem is solved and compared to the Boltzmann solution. Memory ef- 
fects and a slowdown of the relaxation process are the results, they can be partially 
described by using a generalized kinetic equation without quasi-particle approxi- 
mation. 



1. Introduction 

In the past decade the effort invested in ultrarelativistic heavy ion collisions 
(URHIC) has grown considerably H. The general hope is, that at some time in the 
near future one may be able to observe an excursion of strongly interacting matter 
from the state of hadrons before the collision into the phase of a quark-gluon plasma 
(QGP). Consequently, the discussion of possible signals from such a shortlived state 
is quite vivid: Weakly interacting probes like photons or lepton pairs, as well as 
strongly interacting signals like those presented by quark flavors of higher mass 
have been proposed. Similar to most of these investigations is the assumption of a 
thermalized plasma phase, followed by the calculation of the time evolution along 
one or the other line of physical reasoning. 

The present paper is a study of the time scales necessary for such a thermal- 
ization. Ultimately, it is the goal to investigate a physical scenario that one may 
reach in future URHIC: A sea of gluons, initially at low temperature, is heated to a 
very high temperature over a short time. In this hot glue, quark-antiquark pairs are 
popping up - until at the very end a thermal equilibrium in the sense of a degen- 
erate plasma is reached. For the purpose of this conference contribution however, 
the calculations will be presented on a more abstract level. The full consideration 
is the subject of an extended papeiS. 

The primary motivation for this study are serious doubts that the requirements 
for the applicability of standard transport theory ( = kinetic gas theory) are ful- 
filled in a QGP: The thermal scattering of constituents occurs so frequently, that 
subsequent collisions overlap quantum mechanically. This implies, that a treatment 
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in terms of quasi-particles is inadequate, one has to account for a nontrivial spectral 
function of the system components 

The use of a finite temperature field theoretical formulation with continuous 
spectral function is also suggested by the Narnhofer-Thirring theorem^, which states, 
that interacting systems at finite temperature cannot be described by particles with 
a sharp dispersion law. As an additional benefit, this approach is free of unphysical 
infrared singularities occuring in standard perturbation theory. 

The paper is organized as follows: In the next section a brief introduction to the 
formalism necessary for non-equilibrium quantum fields is given. In section 3 the 
approximate spectral function is discussed, followed by a solution of the quantum 
transport equation in section 4. In section 5 a generalized kinetic equation is solved 
which stands between the usual Boltzmann equation and the quantum transport 
equation of section 4. Conclusions are drawn in the final section of the present 
work. 

2. Matrix- valued Schwinger-Dyson equation 

As has been pointed out by various authors, the description of dynamical (time 
dependent) quantum phenomena in a statistical ensemble necessitates a formalism 
with a doubled Hilbert spaceS!. For the present purpose the relevant content of 
this formalism is, that its two-point Green functions are 2x2 matrix- valued. It is 
left to the reader to chose either the conventional Schwinger-Keldysh, or Closed- 
Time Path (CTP) Green function formalism,! or the technically simpler method of 
Thermo Field Dynamics (TFD)0. 

Within this matrix formulation, consider the Schwinger-Dyson equation for the 
full quark propagator S = So + So ©HQS. Here S is the free and S the full two-point 
Green function of the quark field, E is the full self energy and the generalized product 
of these is to be understood as a matrix product (thermal and spinor indices) and an 
integration (each of the matrices is a function of two space coordinates). Throughout 
this paper the convention is used to write space-time and momentum variables also 
as lower indices, e.g. T, xy = Ti(x,y). 

In the CTP formulation as well as in the a = 1 parameterization of TFDEJ, the 
matrix elements of S, S and S obey 

ell i c22 _ Cl2 I c21 y^ll I ^22 _ yil2 yi21 /-i \ 

D (0) + D (0) — D (0) + *->(()) 2^+2^ — —2^ — 2^ . {i) 

Therefore the four components of the Schwinger-Dyson equation are not indepen- 
dent, the matrix equation can be simplified by a linear transformation which one 
may conveniently express as a matrix product It achieves a physical inter- 

pretation only in the TFD formalism, see ref. i. The transformation matrices B 
are 

sw^' 1 -"' -';), (2) 



depending on one parameter only. For example, the third term in the Schwinger 
Dyson equation becomes 



B(n) r s S„ O E S <S(»))- = ( S '' & E " ^ ^ ££f£$ ) . (3) 

Here, r 3 = diag(l, —1), T. R,A are the retarded and advanced full self energy function, 
and S R ' A are the retarded and advanced full propagator (similarly for So) 

S R = s 11 + s 12 , S A = s 11 + s 21 

S R = gll _ S 12 ^ gA = S ll _ S 21 _ ( 4 ) 

The diagonal elements of the transformed equation therefore are retarded and ad- 
vanced Schwinger-Dyson equation. The off-diagonal element is a transport equation. 

Now one switches to the mixed (or^Wigner) representation of functions depend- 
ing on two space-time coordinates: Exp = jd 4 (x — y) exp (iP M (x — yY) T, xy with 
X = (x + y)/2, the -sign will be dropped henceforth. The Wigner transform of the 
convolution S G is a nontrivial step. Formally it may be expressed as a gradient 
expansion 

Jd\x - y) exp (iP„0 - yY) E xz G zy = exp (-iO) t XP G XP . (5) 

O is a 2nd order differential operator acting on both functions appearing behind it, 
OAxpBxp = \ (dxAxpdpBxp — dpAxpdxBxp)- Obviously, this first-order term 
in the application of the infinite-order differential operator exp(— iO) is the Poisson 
bracketu. Henceforth this operator is formally split into cos O— i sin O. Similarly, one 
defines real Dirac matrix-valued functions as real and imaginary part of propagator 
and self energy: 

Sx'p = G X p =F inAxp Exp = Re£ X p =F ivrr^p • (6) 

Axp is the generalized spectral function of the quantum field. 

Now consider the equations obtained by action of Dirac differential operators 
(= inverse free propagators) on the matrix-transformed Schwinger-Dyson equationl! 2 ! 
The diagonal components are 

Tr [(P%, - m) Axp] = cos O Tr [ReS XP A X p + T XP G X p] 

Tr [(P^ t - m) Gxp] = Tr [1] + cos O Tr [ReS XP G X p - vr 2 T XP A XP ] ■ (7) 

Two important facts about these equations have to be emphasized. First notice 
that these equations do not in general admit a 5-function solution for the spectral 
function Axp even in zero order of the gradient expansion. This has led to erro- 
neous statements in papers deriving transport equations from the Schwinger-Dyson 
equationEl, because the right side of (0) may not be disregarded. In short terms, 
there is not such thing as a mass shell constraint in quantum transport theory ! 



Secondly, the equations do not contain odd powers of the differential operator O. 
This implies, that when truncating the Schwinger-Dyson equation to first order in 
this differential operators (the usual order for the approximations leading to kinetic 
equations), the spectral function Axp may still be obtained as the solution of an 
algebraic equation. 

The off-diagonal component of the transformed Schwinger-Dyson equation reads, 
after acting on it with the inverse free propagator a\H 



D D xy 



xz ^ zy 



K 



y K a) ^ A 

^xz w ^ zy i 



with kinetic components S K = (1 — n) S 12 + n S 21 and T, K 
Inserting the real function 

which henceforth is labeled quantum transport equations 



(1-n) S 12 + nE 21 . 

Inserting the real functions defined before, this leads to a differential equation, 

iona^A; 



Tr 



7 M + 2 sin O ReSxp + cos O 2-kT X p) S X p 



2iTr 



i sin O T<xp Gxp — cos O T, xp inAxp 



K 



(9) 



Note, that here even as well as odd powers of the operator O occur. The solution 
in zero order O is not trivial, since it leads to the diagonalization of the propagator 
in equilibrium state*!'". 



3. Effective fermion propagator and spectral functions 

In a thermal equilibrium state at temperature T, the full propagator of a fermionic 
quantum field has to obey the Kubo-Martin-Schwinger boundary conditionBIE: 



1 - n F (E)) S™(E,fi + n F (E)S 2 e l Q (E,p) = . 



eqv 



(10) 



np(E) is the Fermi-Dirac equilibrium distribution function at temperature T, 
np(E) = (e^ B ~ M ) + l) -1 . As seen above, the matrix valued propagator has only 
three independent components, two of which are furthermore complex conjugate. 
One may now use the KMS condition to ehminate the off-diagonal component of 
the equilibrium propagator in favor of n(Ep&3: 



4 a qWp)= JdEA(E,p) x 



r,{B{n{E)))- 1 



Po - E + ie 



V 



Po 



E 



B{n{E)) . 



ie 



Here A(E,p) is the spectral function of the quark field, properly normalized and 
approaching a 5-function for vanishing interaction. 



With the present paper one is addressing non-equilibrium states. For such states 
one may not derive a spectral representation of the propagator in general^, but one 
may still exploit the fact that retarded and advanced propagator are by definition 
analytical functions of the energy parameter in the upper or lower complex energy 
half plane. 

Hence, even for non-equilibrium states one may write in the mixed (or Wigner) 
representation 

oo 

S R ' A (E,p,X) = ReG XP TniAxP= [ dE' A(E' , p, X) —— , (12) 

— oo 

since this is nothing but the Wigner transform of S^ A = =F27ri0 (±(x — yo)) A xy - 

By inspection of eq. (0) one finds, that only a self energy function is needed for a 
full determination of the function Axp- This self energy function is in general a func- 
tional of Axp again - which then leads to a complicated set of integro-differential 
equations for the self consistent determination of the retarded and advanced prop- 
agator. 

For the limited purpose of the present paper however, one makes some physically 
motivated assumptions: 

1. The self energy function for the quarks is dominated by gluonic contribu- 
tions. This is justified because the quark-quark scattering cross section is 
much smaller than the quark-gluon cross section. 

2. The gluon background is dominated by external conditions, i.e., we neglect 
the back-reaction of quarks on the gluon distribution. 

3. The external conditions determining the gluon field are changing in a short 
time interval, and the system is translationally invariant in 3-dimensional 
coordinate space. 

4. One neglects the influence of anti-quarks in the spectral function. This re- 
striction is removed in the extended version of this paper, ref. i. 

For these assumptions also exists a practical reason: They allow a clean separation 
of various aspects of the quantum transport problem, whereas this separation is 
difficult (if not impossible) when considering more realistic systems. 

These assumptions lead to the following ansatz for the imaginary part of the self 
energy function: 

( Ti if -A > t 

ttTxp = 7° T t = 7° gT(t) =7% + A )p ~ ^ jf > t > —A . (13) 

[T f if t > 

Within this ansatz the limit of A — ► is discussed separately, it corresponds to 
instantaneous heating of the gluon background. Furthermore one ensures causality 
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Figure 1: Time dependent spectral width parameter 7 t . 
Parameters are #=0.12, T t =l MeV, T f = 200 MeV, m = 10 MeV. 
Thin lines: T t from eq. (|TB|), thick lines: 7 t from eqs. (15), 
continuous lines: A = 0, dashed lines: A = 4 fm/c. 



by calculating the real part of the self energy through a dispersion integral. This 
integral is divergent, hence in principle one also needs a regularization procedure - 
but the effects of this divergence cancel in the equations. 
For the quark spectral function, one uses the simple form 

A(E,p,t) = ^ - 7 ; - . (i4) 

7T (E - u t ) + Y t 

Hence, one approximates the quark spectral function by two time-dependent param- 
eters u t and 7 i? which may be interpreted as effective mass and effective spectral 
width. One may argue about the validity of this approach, in particular whether 
not a momentum dependent spectral width is an absolute necessity for a realistic 
calculation. 

However, first of all one may safely assume that the quarks appearing in the hot 
medium are slow - hence the properties of the quark distribution may be approxi- 
mated by those of quarks at rest. A second argument in favor of this ansatz is the 
question of causality: The expectation value of the anti- commutator of two quark 
fields is nothing but the Fourier transform of the spectral function. Hence, while 
for some more general spectral function causality may be violated, the above ansatz 
guarantees it when supplemented with a corresponding antiparticle piecei'El. 

With the above spectral function the coupled system ([7]) reduces to a single 
nonlinear equation for j t , plus the condition uj 2 = uj" 2 = p 2 + m 2 . This latter 
condition is more complicated, when the anti-particle piece of the spectral function 



is taken into account!!. The energy parameter is chosen as E = ujq, which yields 
instead of eq. (|7|) as the Schwinger-Dyson equation for the retarded (or advanced) 
two-point function of the quarks: 

lt = gTt + g(T f - T>) Q(t) + g(T f - T t ) ^-±A _ -±-\ e(-t)G(t + A) 

- 9{Tf ~ Tl) (e(t) e- 2 ^ - e(t + A) e- 2 ^ + A )) (15) 

In the limit A — > 0, this becomes even simpler: 

lt = gT t + g(T f - T 4 ) Q(t) (l - e~ 2 ^) (16) 

In Fig. 1, the solution of these equations is plotted in comparison to the time 
dependent imaginary part of the self energy function from eq. (|T3|). It is obvious, 
that the solution of the nonlinear equations (|15|) resp. (|16D approaches the imaginary 
part of the self energy function with a characteristic delay time. Simply using T t 
from eq. (|T3|) instead of 7t - which would correspond to an adiabatic approximation 
- therefore ignores this delay time. In ref.i it is discussed how this delay time is 
calculated from the system parameters. 



4. Transport equation 

As was stated above, the off-diagonal component of the transformed Schwinger- 
Dyson equation is a transport equationOu. To see this more clearly, define the 
generalized covariant distribution function N X p through the equation 

(1 - N XP ) S% P + N XP Sf P = . (17) 

Note the similarity with eq. (0): The above equation indeed ensures, that in the 
limit of thermal equilibrium one achieves lim equ jiAfxp = np{E). For the purpose 
of the present paper Nxp is taken as a scalar function. The description of phe- 
nomena like spin diffusion requires to use a Dirac matrix valued Nxpi^- It follows 
that Sxp = 27ri (Nxp — n) Axp- From this step the mathematical interpretation 
of the generalized distribution function Nxp is obvious: It is the parameter which 
diagonalizes the the full nonr£fluilibrium matrix- valued propagator through the Bo- 
goliubov matrix B from (^EM: 

B(Nxp) r 3 Sxp (B(N xp ))-' = ( ° XP ~ ™ AxP ). (18) 

\ LrXP + WAxP J 

For the following, one furthermore defines a "pseudo-equilibrium" distribution func- 
tion: The 2x2 matrix structure of self energy function allows to diagonalize it also 
by a Bogoliubov transformationta with a parameter N XP such that 

S^ 2 P = 2mN XP T X P E| P = 2vri (n xp - l) T XP . (19) 
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Figure 2: Normalized time dependent fermionic distribution function for slow quarks. 
Parameters as in Fig. 1; thin lines: Nf /np(m,Tf) from the Boltzmann equation 
fl2"Ip, thick lines: N t /n F (m,Tf) from the quantum transport equation (|2"DD; 
continuous lines: A = 0, dashed lines: A = 4 fm/c. 



In the present approach Nxp is determined by the hot gluon gas acting as back- 
ground, hence without the back-reaction it is equal to the equilibrium function, 
N^p = np(E, T(t)) with a time dependence due to the time dependence of the tem- 
perature. Looking at slow quarks with E — m, one furthermore replaces N(X; m,p) 
by Nt and neglects the energy derivative of np{E,T(t)). The resulting quantum 
transport equation according to @ then is: 



with T(t) as defined in eq. ([13|). This equation looks surprisingly similar to a kinetic 
equation in relaxation time approach. However, this similarity is superficial: The 
kinetic equation, or Boltzmann equation, derived for this simple model system reads 



with the imaginary part of the self energy T t from eq. ([13]) instead of the spectral 
width parameter jf. That these differ substantially in the beginning of the relaxation 
process has been shown in the previous section. 

Fig. 2 depicts the influence of this difference on the solution of the transport 
equation. The result is that the relaxation process is slowed by the inclusion of 
the spectral function of the system components. Please observe, that the curves of 
Fig. 2 employ the same behavior as seen in Fig.l: The relaxation rate is similar 




2 lt (N t -n F (m,T(t))) 



(20) 




(21) 



in the quantum transport and the Boltzmann equation, but the former experiences 
a characteristic delay time with respect to the latter. This delay time is almost 
doubled with respect to the delay time occuring in the spectral width parameter j t , 
an asymptotic calculation is carried out in ref.B. 



5. Gradient expansion 

One may now raise the question, whether one can produce an equation which at 
least takes some of the quantum features of particles into account in an otherwise 
kinetic picture. The reason for this is, that in a general non-equilibrium system one 
cannot hope to reduce the equations (0) and @ to such simple forms as obtained 
above. Even a purely numerical solution of these equations seems to be impractical 
if not impossible. 

Therefore, to answer the question, consider the two steps which are between 
the equations ( p0|) and (0]): First of all a quasi-particle approximation, secondly 
an expansion of the operator exp(— iO) to first order, i.e., replacing it by 1 — iO. 
The first of these steps would be in contradiction to the philosophy outlined in the 
introduction to this work. The second step however may be kept: To expand the 
diagonal as well as the off-diagonal pieces of the original matrix- valued Schwinger- 
Dyson equation to first order in the operator OoEI. 

The necessary differential equation for Nxp has been derived in ref.0, correct 
to first order in the gradient expansion it reads 



Tr 



A 



XP 



m 



ReS 



XP 



N 



XP 



iTr 



Axp N X pV$ P - (N XP - 1) 

^XP 



(It 



dE 



sin(rE) Tr 



A{X;P + E,P), 



N XP ^ (t + r/2, X; P) - (N XP - 1) S 12 (t + r/2, X; P) 



N- 



(22) 



In this equation, {-, •} denotes the Poisson bracket, the index N means that the 
derivatives are not acting on N XP . Here, as outlined before, one may use a spectral 
function which is the solution of an algebraic equation. For the present simple 
model this means to replace % by T t in the function A. Note, that the above 
equation is strictly causal: It involves a time integral only over the past history of 
the system, and its derivation is based on the dispersion integral (|i~2|) . The problem 
of unphysical singularities in the propagator therefore does not occur. 

Furthermore, replacing N XP by the unknown function and inserting all the 




Figure 3: Normalized time dependent fermionic distribution function for slow quarks. 
Parameters as in Fig. 1; thin lines: left Nf /n F (m,Tf) from the Boltzmann equa- 
tion (]2"ID, right N t /n F (m,Tf) from the quantum transport equation (pD); thick lines: 
iV t /n F (m,Tf) from the generalized kinetic equation (|23|); continuous lines: A = 0, 
dashed lines: A = 4 fm/c. 



previous definitions, one obtains the nonlinear equation 

jNf = -2T t (N?-n F {m,T{t)) 



+2g{T f -T i , 



W(T + 7^r)exp(-2r t t) 



A 2T t A 



-9(t + A) 



t + A 



+ e(-()e(t+A) 



A 2r t A 

1 



2r,A 



exp(-2r t (t + A)) 
N G _ n F (m,T f )T f - n F (m,Ti)Ti ^ ^ 



In the limit A — > this may be simplified to 

jN? = -2T t (N?- nF (m,T(t)) 

+4t9(t) (giTj-Ti)) 2 exp(-2T t t) 



G n F (m, T f )T f - n F (m, T^)T { 



Tf - 71; 



(24) 



Shown in Fig. 3 is the numerical solution for in comparison to the Boltzmann 
solution Nf as well as the full quantum transport solution N t . 



6. Discussion and Conclusion 



The comparison of the three methods to describe the relaxation problem of a 
quark-gluon plasma (QGP) shows, that the full quantum transport equation results 
in a much slower equilibration process than the Boltzmann equation. This result 
is in agreement with other attempts to solve the quantum relaxation problemiS 
The quantum system exhibits a memory, it behaves in an essentially non-Markovian 
way. 

In particular, for the physical scenario studied here, the system "remembers" 
that it has been equilibrated some time ago. The relaxation rate then is very similar 
to the Boltzmann rate, but the system follows with a characteristic delay time. This 
delay time depends on the system parameters in a non-algebraic way, hence one may 
be subject to surprises for physical examples. 

In the present quantum transport example for the QGP, the time to reach 1- 
1/e 2 ~ 86 % of the equilibrium quark occupation number is almost doubled (14.7 
fm/c as compared to 8.2 fm/c in the Boltzmann case). Thus, it may be carefully 
stated, that the question of the applicability of standard transport theory with 
quasi-particles needs further investigation: It might turn out, that quantum effects 
(= memory as decribed in this contribution) substantially hinder the thermalization 
of a QGP over long time scales. 

One also finds, that this result holds for instantaneous as well as fast (A = 4 
fm/c) heating of the bosonic background. Without elaboration at this point it may 
be stated that the inclusion of antiquarks into the spectral function does not change 
these figures substantially; it only leads to small oscillations of the relaxation rate 
around the value given in Fig. 1. 

The calculated numerical value of 14.7 fm/c for the thermalization time of slow 
quarks is certainly so large, that the cooling of the bosonic background has to be 
taken into account for realistic estimates. Thus however one runs into the principal 
problem of non-equilibrium quantum field theory: The solution of time-dependent 
coupled equations for the Green's functions, hardly possible in any concrete case^ 
A way out of this dilemma might be offered by the generalized kinetic equationMj 
( |2"2] ) , which is related to the quantum transport equation as well as to the Boltzmann 
equation: It does not contain the convolutions over coordinate space that are hidden 
in the Schwinger-Dyson equation. 

However, it does contain the gradient approximation of standard transport the- 
ory - and thus its applicability to the system studied here is questionable, since a 
step function in time certainly involves large gradients. The present comparison is 
therefore justified only through its results: The fact, that with the generalized trans- 
port equation one does at least partially describe the memory effects in a quantum 
system (the characteristic time now is 11.4 fm/c) is encouraging. Applications of 
this transport equation to more complicated systems seem to be possible, at least in 
cases where one previously has used Boltzmann-like or Vlasov-like equations which 
also contain this gradient expansion to first order. 

As a more general remark at the end of this paper it might be added, that the 



present results certainly demonstrate the importance of solving all three components 
of the matrix-valued Schwinger-Dyson equation on the same level of approximation. 
Using only a trivial approximation to the diagonal equations, i.e., replacing the 
spectral functions of the model by some "mass-shell constraint" , is not justified for 
strongly interacting hot systems. 
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